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Jet physics in e + e annihilation from 14 to 209 GeV 

S. Kluth a 

a Max-Planck-Institut fiir Physik, 

Fohringer Ring 6, D-80805 Munich, Germany 

Results from the study of jets in hadronic final states produced in e + e~ annihilation at y/s = 14 to 209 GeV are 
discussed. Precise measurements of the strong coupling as at many points of y^s provide convincing evidence for 
the running of as as expected by QCD. Several measurements of the running b-quark mass mb(M z o) are combined 
yielding mb(A/ z o) = (2.92±0.03(stat.) ±0.31(syst.)) GeV and compared with low energy results. Strong evidence 
for the running of the b-quark mass is found. Experimental investigations of the gauge structure of QCD are 
reviewed and then combined with the results for the colour factors C a = 2.89 ± 0.03(stat.) ± 0.21(syst.) and 
C ' f = 1-30 ± 0.01(stat.) ± 0.09(syst.) and correlation coefficient p — 0.82. The uncertainties on the colour factors 
are below 10% and the agreement with the QCD expectation Ca = 3 and Cf = 4/3 is excellent. 



1. INTRODUCTION 

Hadron production in e + e _ annihilation is an 
ideal laboratory for QCD studies. For hadronic 
final states there is no interference with the lep- 
tonic intial state. The energy of the hadronic final 
state is maximal in the laboratory system (in the 
case of symmetric colliders) thus allowing efficient 
production and experimental study of highly en- 
ergetic hadronic systems. Since the particles in- 
volved in the electroweak production of hadrons 
via the process e + c~ — * (Z°7)* — > qq are point- 
like there are no parton density functions to take 
into account. 

The current understanding of hadron produc- 
tion in e + e~ annihilation in the Standard Model 
is shown schematically in figure ^ The collid- 
ing electron and positron annihilate into a pho- 
ton or Z°, possibly with additional radiation of 
energetic photons from the incoming beam par- 
ticles (initial state radiation or ISR). The inter- 
mediate photon is always virtual with its invari- 



ant mass equal to the invariant mass v s' of the 
e + e~ system at collision after ISR while the Z° 



can be produced on mass-shell when ys'~ M z o . 
The intermediate photon or Z° then decays into 
a qq state which can in turn radiate gluons. The 
continuing processes of gluon radiation and qq 
production result in a parton shower which only 



Electro-weak Production Parton Shower Hadronisation 




d 1 

u J 

and many more 
u 

d } A°^Jl p + 

S 

d _ } d^kV 



Vs' - 1 GeV 

Parton Level Hadron Level 



Figure 1 . Schematic view of hadron production 
in e + e _ annihilation. 



terminates when the partons (quarks and gluons) 
reach energy scales of ~ 1 GeV. Final states at 
this stage are referred to as parton-level. At these 
scales the increasing strong coupling causes the 
formation of colourless bound states of quarks, 
i.e. hadrons. This process, known as hadroni- 
sation, is not well understood at a fundamental 
level in QCD but many successful models based 
on QCD have been developed. The hadronic final 
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state after decays of short lived particles (usually 
denned to have lifetimes < 300 ps) is said as being 
at the hadron-level. 

It is well known that the hadronic final states in 
e + e _ annihilation show a so-called jet structure, 
i.e. several groups of particles are observed in 
relatively small and separate regions of the solid 
angle around the interaction point. These struc- 
tures are interpreted as remnants of the qq pair 
from the virtual photon or Z° decay and possi- 
ble energetic gluon radiation. Many definitions 
of jet finding algorithms |1I2| or event shape ob- 
servables (see e.g. |3I4| 1 are used to quantify the 
jet structure and compare with predictions from 
QCD. 

Jet finding algorithms cluster particles lying 
closely together in phase space. As an exam- 
ple in the Durham algorithm [S] the distance 
between two particles i and j with energies E^ 
and Ej and spanning the angle 9ij is defined as 
Uij = 2EiEj(l — cos 0ij)/s. The pair of particles 
with the smallest i/ij is combined into a pseudo- 
particle by adding the four-momenta and is then 
discarded. The procedure is iterated until only 
one jet remains. One can now study the jet struc- 
ture as a function of the resolution parameter y cut 
and e.g. determine the fraction of 3-jet events 
at a fixed value of y cu t- Alternatively, the value 
of yij = 1J32 where an event changes from a 3- 
jet to a 2-jet configuration is characteristic for its 
structure, i.e. a large value of z/32 indicates the 
presence of a significant third jet. 

At the time of writing data of observables such 
as jet production rates or event shape distri- 
butions measured at centre-of-mass energies y/s 
ranging from 14 to 209 GeV are available. These 
are especially interesting, because many QCD 
predictions are strongly energy dependent due to 
the running of the strong coupling ag(Q). In 
leading order the running of as is expressed as 

as(Q) = as(fi)/(l + 2f3 a s (ji) ln(j*/Q)) (1) 

with p Q = (IXC A - 2n/)/(127r). The value a s (Q) 
at a scale Q is given by as (p) at a reference scale 
\x and the logarithmic ratio of the energy scales 
21n(Q//x). However, hadronisation effects scale 
like 1/Q n , with n = 1 for many jet algorithms or 
event shape observables, implying that a mean- 



ingful comparison of QCD predictions with data 
requires a detailed understanding of the scale de- 
pendencies of perturbative QCD and as well as of 
the hadronisation effects. The present e + e~ data 
allow to study these problems in a controlled way. 
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Figure 2. Cross sections for various e + e anni- 
hilation processes as indicated on the figure UJ. 



The selection of c + e~ annihilation into 
hadronic final states has to consider different ex- 
perimental conditions at the various cms ener- 
gies. Figure |2] shows the total cross sections for 
hadron production, pair production and 

two-photon interactions [H]- At y/s < M^o back- 
grounds stem from two-photon processes, t + t~ 
pair production and events with hard ISR. For 
example in the JADE data recorded between 
y/s = 14 and 44 GeV after demanding large to- 
tal energy, balance of momentum along the beam 
axis and large particle multiplicity the total back- 
ground is reduced to about 1% with a selection 
efficiency of about 60% 7 . On the Z° peak simi- 
lar event selections reduce the background to less 
than 1% with almost 100% selection efficiency. 
Above the Z° peak for y/s > 2Af\y backgrounds 
from hadron production mainly via production of 
W + W~ pairs decaying both hadronically become 
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important. The prominent 4-jet structure of the 
W + W~ events and their different kinematic prop- 
erties compared to e + e~ — > qq events allows to 
reduce backgrounds to less than 10% at selection 
efficiencies of about 80%. 

2. HARD QCD 

2.1. Measurements of as 

The strong coupling is in QCD predictions for 
massless partons the only free parameter. A fit of 
the theoretical prediction to the data with as &s 
a free parameter yields a measurement of ag(Q) 
where Q is identified with the hard scale of the 
process, i.e. Q = y/s. Hadronisation effects are 
not part of the perturbative QCD prediction and 
are usually treated as a correction which is ap- 
plied to the theory before comparison with the 
data. The choice between several hadronisation 
models implemented in Monte Carlo generator 
programs and uncertainties in the experimentally 
adjusted free parameters of the models are re- 
flected in the systematic uncertainties of the as 
determinations. 

Figure [21 shows the measurement of 
as (206.2 GeV) by L3 using a set of five event 
shape observables Thrust T, Heavy Jet Mass p, 
Total and Wide Jet Broadening Bt and %, and 
C-parameter C. For these observables matched 
C(a^)+NLLA QCD calculations |9ll0lll| are 
available, which provide reliable predictions 
over a large part of the phase space populated 
by 2- and 3-jet events. The combined result 
of the five measurements is as (206.2 GeV) = 
O.f 13 ± 0.002(stat.) ± 0.005(syst.). 

Figure 0] presents the results of similar anal- 
ysis performed by L3 at the other LEP 2 en- 
ergy points, at LEP 1 on the Z° peak and us- 
ing hadronic Z° decays with hard final state 
photon radiation (FSR) for measurements at 
\fs < Mz'>. The error bars show the exper- 
imental and statistical uncertainties only. The 
solid line shows a fit of the QCD prediction for 
the running of as with the result ag(Mzo) = 
0.123 ± 0.001 exp±0.006(theo.) and x 2 /d-o.f. = 
18/15. The assumption of as=const. results in 
X 2 /d.o.f. = 52/15; the data thus support the run- 
ning of as as predicted by QCD. Similar analyses 
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Figure 3. Data from L3 of event shape observ- 
ables measured at = 206.2 GeV. Superim- 
posed are C(a^)+NLLA QCD fits 0. 



by the LEP collaborations are |12I13| . 

The data of the former JADE experiment taken 
between 1979 and 1986 at PETRA/DESY at 
y/s — 14 to 44 GeV have been re-analysed re- 
cently 7 . The analysis employs modern Monte 
Carlo event generators together with the original 
detector simulation program to obtain more re- 
liable experimental corrections with reduced sta- 
tistical and systematic uncertainties. The contri- 
bution of e + e~ — > bb events is subtracted to re- 
duce the effects of heavy quark masses and heavy 
hadron decays on the analysis. The data can 
be compared directly with QCD predictions for 
massless partons corrected for hadronisation ef- 
fects. Several event shape observables are used to 
determine as in the same way as the LEP collabo- 
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Figure 4. Results of ag(\/s) from L3 compared 
with the expectation from QCD for the running 
of as (solid line) and with as = const, (dashed 
line) [Hj. 



rations using the best currently available matched 
C(a|)+NLLA QCD calculations. 

The results of fits to 1 — T are shown in fig- 
ure [5] The diagrams show the data for 1 — T 
at *fs = 14, 22, 35 and 44 GeV corrected for 
experimental effects. A decrease of the shoulder 
at large 1 — T with increasing y/s is clearly visi- 
ble; this observation corresponds to the expected 
lower fraction of events with hard gluon radiation 
due to the running of as- Superimposed are the 
results of QCD fits using matched 0(a|)+NLLA 
(ln(R)- and R-matching) or C(a§) calculations 
corrected for hadronisation effects, which gener- 
ally describe the data well. The hadronisation 
corrections Ki with uncertainties are shown as 
grey bands. The hadronisation corrections are 
large with large uncertainties at low y/s but be- 
come much smaller with increasing \fs. The re- 
sults for a s (y/s) are a s (14 GeV) = 0.170 ±0.019, 
a s (22 GeV) = 0.151 ± 0.013, a s (35 GeV) = 
0.146 ± 0.012 and a s (44 GeV) = 0.131 ± 0.009. 
The data together with other compareable results 
from the LEP experiments and SLD are in good 
agreement within the experimental errors with 
the QCD prediction for the running as based on 
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Figure 5. Data for 1 - T from JADE at */s = 14, 
22, 35 and 44 GeV are shown on the left with 
C(a^)+NLLA or 0{aQ QCD fits superimposed 
(lines). On the right the hadronisation correc- 
tions K% are shown with uncertainties as grey 
bands 0. 



a s (M z o) = 0.118 ± 0.003 QI|. As above a fit as- 
suming as=const. results in a \ 2 probability of 
- 10" 5 . 

A novel method to determine the strong cou- 
pling as uses the rate of 4-jet events R4 de- 
termined using the Durham algorithm |15I16| . 
Recently the NLO corrections for observables 
sensitive to at least four partons in the final 
state, i.e. corrections of C(a|), became avail- 
able |17I18| . These calculations are matched with 
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existing NLLA calculations. 

A measurement of as (M z o ) based on these pre- 
dictions is presented in figure . The un- 
corrected data for R4 measured using 2.5 • 10 6 
hadronic Z° decays using the Durham algorithm 
are shown as a function of the jet resolution 
parameter y C ut- The line shows the fit of the 
0(a|)-|-NLLA (R-matching) theory corrected for 
hadronisation and detector effects. The fit is per- 
formed using only the bins within the fit range as 
indicated where the total correction factors devi- 
ate from unity by less than 10%. The result of the 
fit for fixed rcnormalisation scale parameter x M = 
1 is a s (M z o) = 0.1170± 0.0022 using the conser- 
vatively estimated total error from ^H] . The fit is 
seen to describe the precise data fairly well. The 
similar analysis of (T^J uses only the 0(a|) calcu- 
lations with experimentally optimised renormali- 
sation scale and finds a s (M z o) = 0.1178 ± 0.0029 
with x 2 , — 0.015. These measurements have com- 
paratively small uncertainties and belong to the 
group of the most precise determinations of as- 

2.2. Running b quark mass 

The QCD predictions discussed so far have 
been made for massless partons. The effects of 
massive quarks on the perturbative QCD predic- 
tions are known in C(a|) as well |19I20I21| . The 
mass of the heavy quark identified with the b 
quark at LEP energies becomes an additional pa- 
rameter of the theory and is subject to renormal- 
isation analogously to the strong coupling con- 
stant. As a result the heavy quark mass nib will 
depend on the scale of the hard process in which 
the heavy quark participates. In leading order 
in the MS renormalisation scheme the running 
b-quark mass is m h (Q) — M h (as{Q)/n) 12 / 23 , 
where Mb is the so-called pole mass defined by 
the pole of the renormalised heavy quark prop- 
agator. The running is known to four-loop 
accuracy |22| . Using low energy measurements of 
TUbi'Tih) one finds that TOb(Mzo) ~ 3 GeV. 

The mass of the b quark is experimentally 
accessible in hadronic Z° decays, because i) about 
21% of the decays are Z° — > bb and these events 
can be identified with high efficiency and purity, 
ii) event samples of O(10 6 ) events are available 
and iii) observables like the 3-jet rate R 3 (y C ut) 
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Figure 6. The 4-jet fraction R4 determined using 
the Durham algorithm is shown as a function of 
J/cut without corrections. Superimposed is a fit of 
C(a|)+NLLA QCD (solid line) [TSj. 



in the JADE or Durham algorithm can have en- 
hanced mass effects going like m 2 /M^ /y cut |23) . 
In the experimental analyses |24I25I26I27I28| the 
ratio B 3 = R b 3 /R 1 3 , R b 3 (R 1 3 ) are the 3-jet rates 
in b quark (light quark) events, is studied which 
reduces the influence of common systematics. 

Figure [7| shows B 3 (y cut ) measured by 
OPAL 1231 before and after detector and hadro- 
nisation corrections for the JADE E0 and the 
Durham algorithm. One observes that the effect 
of the b quark mass is about 5% on B 3 with rea- 
sonable uncertainties as indicated by the bands. 
The effect goes in opposite directions for JADE 
E0 and the Durham algorithm, because in the 
JADE type algorithms based on invariant mass 
to cluster jets the large b quark mass causes an 
enhancement of B 3 while in the Durham algo- 
rithm the reduced phase space for hard gluon 
emmission dominates thus reducing B 3 |26| . 

Figure [8] gives an overview over presently avail- 
able measurements of mb(Mzo) |28|. These 
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Figure 7. The figures show distributions of B3 
determined using the Durham or JADE EO jet 
finding algorithms before and after all corrections. 
The solid lines in the lower row of figures give the 
uncertainties on B3. The vertical dashed lines 
indicate the values of y cu t where the QCD pre- 
dictions are calculated |2*5] . 



are compared with low energy measurements 
and the QCD prediction for the running of 
frit, (Q). The predicted range of mb(M z o) is in 
good agreement with the measurements. The 
measurements of m\ ) (M z o) of 25 26 27 28 can 
be averaged assuming the statistical errors to 
be uncorrelated, the experimental and hadro- 
nisation uncertainties to be partially correlated 
and the theory uncertainties to be fully cor- 
related 1 . The result is m h (M z o) = (2.92 ± 
0.03(stat.)±0.31(syst.)) GeV, which may be com- 
pared with the average of low energy measure- 
ments TO b (m b ) = (4.24 ±0.11) GeV (2H|. The 
difference becomes m\ ) (M z o) — 171^(111^) = (1.32 ± 
0.33) GeV which is non-zero by four standard de- 

1 For partial correlation the covariance matrix element is 
Vx,i = min(o-i , o-j ) 
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defined as Vi.j = min(<r; , rxj ) 2 and for full correlation as 



Figure 8. Measurements of rrib(Q) at low and 
high energies Q are compared with the QCD ex- 
pectation for the running of m,b(Q) |28|. 



viations thus providing strong experimental evi- 
dence for the running heavy quark mass in QCD. 

3. QCD GAUGE STRUCTURE 

QCD as the gauge theory of strong interactions 
assumes that quarks carry one out of three strong 
charges, referred to as colour. The requirement 
of local gauge symmetry under SU(3) transfor- 
mations in the colour space generates the gauge 
bosons of QCD: an oktett of gluons each carrying 
colour charge and anti-charge, see e.g. [23j- The 
gluons can thus interact with themselves; it is due 
to this property of the QCD gauge bosons that 
the theory explains confinement and asymptotic 
freedom via the running of the strong coupling 
a S - 

In perturbative QCD at NLO three fundamen- 
tal vertices contribute: i) the quark-gluon vertex 
with colour factor CU, ii) the gluon-gluon ver- 
tex with colour factor Ca and iii) the qq produc- 
tion from a gluon with colour factor TpUf, see 
e.g. The colour factors are Cf — 4/3, Ca — 
3 and T F n f = 1/2 • 5 in QCD with SU(3) gauge 



7 



symmetry and specify the relative contribution 
of the corresponding vertex to observables. In 
NLO QCD the prediction for an observable R is 
R = Aa s + (B Cf C f + B Ca C a + B TF T F n f ) C F ; 
NLLA predictions e.g. for event shapes or jet 
rates decompose in a similar way. For NLO pre- 
dictions for 4-jet observables an analogous decom- 
position in terms of the six possible products of 
two out of the three colour factors holds |32| . 

Experimental investigations of the gauge struc- 
ture of QCD are possible because of the different 
angular momenta in the initial and final states of 
the fundamental vertices. It is an important test 
of QCD to probe the gauge structure in exper- 
iments. Several techniques with rather different 
experimental and theoretical uncertainties have 
been developed; we will discuss here some recent 
results. 

3.1. Four-jet events 

The LEP experiments ALEPH Q2] and 
OPAL have analysed 4-jet final states from 
hadronic Z° decays using the recent QCD NLO 
predictions (see e.g. and references therein). 
The 4-jet final states are selected by cluster- 
ing events using the Durham algorithm 5 with 
Vcut = 0.008 and demanding four jets. At this 
value of y cu t the 4-jet fraction is relatively large 
(R 4 ~ 7%) and the four jets are well separated. 

The energy-ordered 4-momenta Pi, i = 1, . . . , 4 
of the jets are used to calculate the angular corre- 
lation observables (see |15I33| for details). As an 
example, the Bengtsson-Zerwas angle xbz is de- 
fined by xbz = Z([pi xp 2 ], [ph xpi]), i.e. the angle 
between the two planes spanned by the momen- 
tum vector pairs (pi,P2) and (pa, pi)- Assuming 
that energy ordering selected the primary quarks 
from the Z° decay as (j>i,P2) the observable %bz 
is sensitive to the decay of an intermediate gluon 
to gluons (vertex ii)) or quarks (vertex iii)) and 
the competing process of radiation of a second 
gluon from a primary quark (vertex i)). 

Figure presents the uncorrected distribution 
of cos(xbz) measured by ALEPH JS]- Superim- 
posed on the data points is the result of a simul- 
taneous fit of the NLO QCD predictions to four 
angular correlations including xbz and the distri- 
bution of R4. From the fit values for as (Mz° ) and 
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Figure 9. The distribution of |cosxbz| before 
corrections is shown (solid points). Superimposed 
is a fit of NLO QCD (solid lines) PS] . 



the colour factors Ca and C F are extracted. The 
results from ALEPH are a s (M z o) = 0.119±0.027, 
C A = 2.93 ± 0.60 and C F = 1.35 ± 0.27; from 
OPAL we have a s (M z o) = 0.120 ± 0.023, C A = 
3.02 ± 0.55 and C F = 1.34 ± 0.26. The system- 
atic errors are dominated by uncertainties from 
the hadronisation corrections and theoretical er- 
rors from estimating missing higher order contri- 
butions 2 . 

3.2. Event shape fits 

In this analysis the decomposition of the 
C(ag)+NLLA QCD predictions for event shape 
observables into terms proportional to the colour 
factors is used. Since the sensitivity of event 
shape distributions measured at LEP 1 alone is 
not sufficient data from y/s — 14 to 189 GeV 
are used. In this way the colour structure of 

2 Both analyses use an unconventional method of estimat- 
ing systematic uncertainties. A more conservative evalu- 
ation leads to total errors for ag(M z o), Ca and Cp of 
±0.048, ±1.06 and ±0.46 for ALEPH and ±0.049, ±1.07 
and ±0.47 for OPAL. 
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the running of the strong coupling contributes as 
well. Hadronisation corrections are implemented 
using analytic predictions known as power correc- 
tions (see US] and references therein). The power 
corrections essentially predict that hadronisation 
effects cause a shift of the perturbative distri- 
bution proportial to 1/y/s in the region where 
the C(ag)+NLLA predictions are reliable. The 
advantage of using power corrections instead of 
Monte Carlo model based hadronisation correc- 
tions is that the colour structure of the power 
corrections is known and can be varied in the fit. 

In the analysis simultaneous fits of as(Mzo), 
Ca and Cf to data for the event shape observ- 
ables 1 - T at Vs = 14 to 189 GeV and C 
at y/s — 35 to 189 GeV are performed. The 
data for 1 — T and C are analysed separatly 
and the results are combined. The results are 
a s (Af z o) = 0.119 ± 0.010, C A = 2.84 ± 0.24 and 
Cf = 1.29 ± 0.18 and are shown on figure IT2lbe- 
low. The errors are dominated by uncertainties 
from the hadronisation correction and from ex- 
perimental effects. 

3.3. Scaling violation in gluon and quark 
jets 

Jets originating from quarks or gluons should 
have different properties due to the different 
colour charge carried by quarks or gluons |36I37| . 
In an analysis by DELPHI [jJH| the scaling viola- 
tion of the fragmentation function (FF) in gluon 
and quark jets at different energies is compared. 
From a sample of 3.7 • 10 6 hadronic Z° decays 
planar 3-jet events with well reconstructed jets 
are selected using the Durham or Cambridge al- 
gorithms, but without imposing a fixed value of 
?/ cu t in order to reduce biases In events 

with the two angles between the most energetic 
jet and the other jets in the range between 100° 
and 170° a b-tagging procedure is applied to the 
jets. Gluon jets are identified indirectly as the 
jets without a successful b-tag. Jets originating 
from light (udsc) quarks are taken from events 
which failed the b-tagging. A correction proce- 
dure based on efficiencies and purities determined 
by Monte Carlo simulation yields results for pure 
gluon or quark jets. 

The jet energy scale is calculated according to 



«h = Ej Ct sin(0/2), where 9 is the angle w.r.t. 
to the closest jet. This definition takes colour 
coherence effects into account 00] • The FF in 
a jet D H (xe, kh) is given by the distribution of 
xe = -Ehadron/ for the hadrons assigned to the 
jet. The evolution of the FFs of gluon and quark 
jets with jet energy scale kh is studied in in- 
tervals of Xb- Figure 1101 presents the quantity 
dD H (xe, KH)/dlog(KH), i.e. the slopes of the 
scaling violation for a given interval of x-e, for 
quark and gluon jets. The steeper slopes corre- 
sponding to stronger scaling violations of gluon 
compared to quark jets is clearly visible. A fit of 
the scaling violations based on the LO QCD pre- 
diction (DGLAP equation) allows to extract the 
ratio Ca/Cf resulting in C A /C F = 2.26 ± 0.16. 
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Figure 10. Measurements of 

dD H (xe, KH)/dlog(Kn) are shown for gluon 
(solid points) and quark jets (open circles and 
solid triangles). The lines represent a fit of the 
LO QCD prediction gHj- 



3.4. Multiplicity in unbiased gluon jets 

In this analysis 0T] the evolution with energy 
scale of the charged particle multiplicities 
and for gluon and quark jets is used to mea- 
sure the colour charge ratio Ca/Cf- The inclu- 
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sive multiplicity iV-Sj is measured in hadronic Z° 
decays as a function of the 3-jet transverse mo- 
mentum scale fcx,Lu — m ( s qg s qg/( s ^QCD) where 
s qg and Sq g are the invariant masses between the 
gluon and the quark jets |39I42| . The analysis 
reconstructs three jets in all events by adjusting 
the value of y cut in the Durham algorithm, se- 
lects one-fold symmetric Y-events |43| and iden- 
tifies the lowest energy jet as the gluon jet. By 
comparing with measurements of N^(k±^ u ) the 
multiplicity in unbiased gluon jets N^' can be ex- 
tracted as a function of the jet energy scale k± Lu'- 
N^-Qcxm) = 2(^ h cl .(fc±, Lu ) ~ Ng(k ±M )). 

Figure^Jshows the data for iV~' together with 
compareable data at lower and higher jet energy 
scales fcj^Lu from CLEO and OPAL. The solid 
line represents the result of a fit of the QCD 
prediction for the evolution of N^' compared to 
|HJ. The fit with Ca/C f as a free param- 
eter yields C A /C F = 2.23 ± 0.14. This analysis 
is complementary to the analysis of scaling viola- 
tions discussed above because the multiplicity is 
dominated by low energy particles while the scal- 
ing violation is dominated by high energy parti- 
cles with large ice- 




Figure 11. Measurements of N^-(Q) are shown 
using 3-jet events (solid points) and other data 
(solid triangles). Superimposed are predictions of 
Monte Carlo generators (dashed and dash-dotted 
lines) and of the QCD fit (solid line) gT| 



3.5. Colour factor averages 

The measurements of the colour factors Ca and 
Cf or of Ca/Cf — x discussed above can be com- 
bined into average values of Ca and C F taking 
into account correlations between Ca and Cf as 
well as between different experiments. The vari- 
ables x and y — Tf/Cf are used to define the % 2 
function 

x 2 = Y.^-^ y] ^-y)+ ( 2 ) 

i 

^2(xi - x)v^f (Xj -x) + 

^2(yi - y)4f (vj -v) > 

where x and y are the averages, the indices i,j 
count experiments, v\ xy \ and are ele- 
ments of the inverses of the corresponding covari- 
ance matrices for the Xi and yi within experiment 
i and for the Xi and yi between experiments. The 
averages x and y are converted to the average 
values Ca and C f after the fit is performed. 

The input data for x and y are directly taken 
from 15 33 38 41 while the results from have 
to be converted 3 . 

The covariance matrices between experiments 
are constructed as follows: ALEPH and OPAL 
4-jet analyses experimental and hadronisation er- 
rors partially and theory errors fully correlated, 
4-jet analyses and event shape analysis hadroni- 
sation errors partially and theory errors fully cor- 
related, and 4-jet and event shape analyses and 
DELPHI FF and OPAL analyses theory er- 
rors partially correlated. 

The averaging fit is done using only the first 
term of equation[21to determine the averages Ca 
and Cf, using only statistical correlations in the 
first term to determine the statistical errors and 
using the full covariance matrix to determine the 
total errors. The final results are Ca — 2.89 ± 
0.03(stat.)±0.21(syst.), C F = 1.30±0.01(stat.)± 
0.09(syst.) with correlation coefficient p = 0.82. 
The relative total uncertainties are about 8% for 
both Ca and C p- 

3 The results are x = 2.20 ± 0.21(stat.) ± 0.25(syst.), 
y = 0.388 ± 0.021(stat.) ± 0.051(syst.), p stat . = 0.98 and 

Psyst. = 0.86. 
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Figure El presents the results of the individ- 
ual analysis in a CV vs. Ca plane together with 
the combined result and the expectations of QCD 
based on the SU(3) gauge symmetry and various 
other gauge symmetries. The correlation coeffi- 
cients for |15I33| were calculated from the refer- 
ences 4 . The error ellipses refer to 86% CL. The 
combined result is in good agreement with the in- 
dividual analyses and with standard SU(3) QCD 
while the total uncertainties are substantially re- 
duced. The other possibilities for gauge symme- 
tries shown on the figure are clearly ruled out. 

4. CONCLUSIONS 

Jet physics in e + e~ annihilation based on data 
from y/s = 14 to 207 GeV allows to study many 
aspects of QCD. Precision measurements of the 
strong coupling as at many points of y/s gave 
convincing evidence for the running of the strong 
coupling as predicted by the theory. 

A combination of measurements of the mass of 
the b-quark mb(-Mzo) at the M^o energy scale 
using jet production rates was performed. The 
resulting value 

m b (M z o) = (2.92±0.03(stat.)±0.31(syst.))GeV(3) 

compared with low energy measurements resulted 
in strong evidence for the running of the b-quark 
mass analogously to the running of as ■ 

The investigation of the gauge structure of 
QCD was discussed for several different meth- 
ods: angular correlations in 4-jet final states from 
hadronic Z° decays, global fits of event shape data 
at many points of y/s, the scaling violation of the 
FF of gluon and quark jets and the evolution with 
energy scale of the charged particle multiplicity 
jV|g- determined from 3-jet events. The results of 
the analyses were combined taking correlations 
between the colour factor measurements into ac- 
count with the results: 

C A = 2.89±0.03(stat.)±0.21(syst.) , (4) 
C F = 1.30±0.01(stat.)±0.09(syst.) , 
p = 0.82 . 



4 ALEPH: p = 0.97, OPAL: p = 0.93. 




Figure 12. The figure presents various measure- 
ments of the colour factors Ca and Cf discussed 
in this report. The ellipses show the correlated 
measurements using 4-jet events |15I33| or event 
shape distributions [32] while the lines repre- 
sent the results of determinations of Ca/Cf from 
DELPHI |2H] (dashed) and OPAL EH] (solid). 
The upper solid and dashed lines overlap. The 
grey filled ellipsis displays the combined result 
for Ca and Cp (see text). The solid triangle 
and squares show the expectations for various as- 
sumptions for the gauge symmetry of QCD as 
indicated on the figure. 



The combined results are in good agreement with 
the individual measurements and have substan- 
tially reduced total errors of less than 10% for 
both C a and Cf- The measurements are also in 
good agreement with the expectation from QCD 
C A = 3 and C F = 4/3. 

The study of jets in QCD with e + e~ annihila- 
tion is in very good shape and will continue to 
provide interesting and important new results. 
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